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I.  INTRODUCTION 


During  the  first  year  of  this  grant,  4  papers  have  been  completed  or  are  nearing 
completion  based  on  full  or  partial  support  of  AFOSR  49620-92-J-0498,  including  the 
work  of  two  graduate  students,  as  follows. 


Krishnamurti,  R.,  1993:  Generation  of  Low  Frequency  Modes  in  Turbulent  Rayleigh- 
Bemard  Convection.  To  be  submitted  to  J.  Fluid  Mech. 

Li,  L.  and  T.  R.  Nathan,  1993:  The  global  atmospheric  response  to  localized  time  varying 
forcing:  Zonally  averaged  basic  states.  To  be  submitted  to  J.  Atmos.  Sci. 

Nathan,  T.  R.  and  A.  Barcilon,  1993:  Low  frequency  oscillations  of  forced  barotropic  flow. 
To  be  published  in  J.  Atmos.  Sci. 

Zhu,  K.,  I.  M.  Navon  and  X.  Zou,  1993:  Variational  Data  Assimilation  Experiments  with 
Variable  Resolution  Finite- Element  Shallow- Water  Equations  Model.  Submitted  to 
Mon.  Wea.  Rev. 


Section  II  contains  abstracts  of  these  papers.  Section  III  contains  additional  information 
concerning  some  of  the  areas  of  accomplishment. 
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II.  ABSTRACTS  OF  PAPERS 


Generation  of  Low  Frequency  Modes  in  Turbulent 
Rayleigh-Bernard  Convection 

by 

Ruby  Krishnamurti 
Geophysical  Fluid  Dynamics  Institute 
and  Department  of  Oceanography 
Florida  State  University 
Tallahassee  FL  32306-3017 

Abstract 

Laboratory  experiments  on  turbulent  Rayleigh-Bernard  convection  have  shown  large 
scale  organization  of  two  kinds,  the  first  a  steady  large  scale  flow,  the  second  locally 
periodic  travelling  clusters  of  plumes.  For  Prandtl  number  Pr  =  7,  the  first  of  these  sets 
in  at  Rayleigh  number  R  =  2  x  106  with  the  flow  in  one  direction  along  the  bottom  and  in 
the  opposite  direction  along  the  top  of  the  layer  (Krishnamurti  and  Howard  1981).  Unlike 
the  cellular  flow  at  low  R  which  has  horizontal  scale  comparable  to  the  layer  depth  this 
flow  fills  the  entire  width  of  the  tank  (which  may  be  50  times  the  depth),  or  may  circulate 
around  an  annular  tank.  All  the  plumes  become  organized  with  the  same  tilt,  the  hot 
plumes  from  the  bottom  boundary  layer  drifting  in  one  direction  with  the  large  scale  flow, 
the  cold  plumes  from  the  top  drifting  in  the  opposite  direction.  A  measurement  of  the 
turbulent  velocities  and  their  correlations  showed  that  (Krishnamurti  and  Howard  1983) 
the  Reynolds  stress  divergence  balances  the  viscous  force  on  the  large  scale  (horizontally 
averaged)  flow  u: 

d _ 

dzuw~vdz * 

The  tilt  of  the  plumes  is  such  that  the  momentum  transport  by  the  Reynolds  stress  uw 
is  counter  to  the  transport  by  the  viscous  stress  on  the  mean  flow.  In  a  low  order  mathe¬ 
matical  model  (Howard  and  Krishnamurti  1986)  it  was  shown  that  tilted  flows  can  result 
from  a  symmetry  breaking  bifurcation  from  upright  cellular  flow  to  one  that  tilts  to  the 
right  or  to  the  left. 

Even  more  remarkable  is  the  second  kind  of  organized  flow  in  which  this  large  scale 
flow  becomes  locally  time  periodic  for  R  >  107.  There  is  a  spontaneous  generation  of  low 
frequency  modes.  These  low  frequency  modes  are  found  both  in  Hovmuller  diagrams  and 
in  power  spectra  of  temperature  time  series.  The  Hovmuller  diagrams  shows  that  clusters 
of  tilted  plumes  travel  as  a  coherent  packet  across  the  convecting  layer.  These  packets  of 
plumes  are  followed  by  quiescent  regions  which  are  followed  by  yet  another  packet  and  so 
on.  A  mathematical  model  will  be  presented  in  which  passage  of  a  plume  has  depleted  the 
thermal  boundary  layer  in  a  history-dependent  way:  the  boundary  layer  is  very  thin  where 
the  plume  has  just  passed,  but  has  had  time  to  thicken  where  the  plume  passed  some  time 
ago.  Thus  at  any  moment  there  is  a  thermal  boundary  layer  of  varying  thickness,  and  the 
pressure  gradient  in  it  drives  a  flow  towards  the  thicker  regions.  Taking  this  to  be  slow 
viscous  flow,  there  is  then  a  shear  and  a  down-gradient  momentum  flux  at  the  wall.  This 
boundary  layer  is  supposed  to  erupt  when  and  where  its  thickness  reaches  some  critical 
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value.  Then  assuming  that  the  interior  momentum  flux  is  entirely  by  Reynolds  stresses, 
we  show  by  matching  interior  and  boundary  layer  heat  and  momentum  fluxes,  that  plumes 
move  away  from  regions  of  thick  boundary  layer  at  just  such  a  rate  that  the  ever- thickening 
boundary  layer  keeps  the  plume  “fed”  with  buoyant  fluid.  The  predicted  period  will  be 
shown  to  be  in  reasonable  agreement  with  that  obtained  in  the  laboratory  experiments. 
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Low  Frequency  Oscillations  of  Forced  Barotropic  Flow 

by 

Terrence  R.  Nathan 
Atmospheric  Science  Program 
University  of  California 
Davis,  CA  95616 

and 

Albert  Barcilon 

Department  of  Meteorology  and  Geophysical  Fluid  Dynamics  Institute 

Florida  State  University 
Tallahassee,  FL  32306-3017 

Submitted  to 

Journal  of  the  Atmospheric  Sciences 
Notes  and  Correspondence 
August  21,  1992 
Revised 
July  8,  1993 

Abstract 

Jin  and  Ghil  (1990)  demonstrate  that  for  topographically  resonant  flow,  low  frequency 
finite-amplitude  oscillations  may  arise  from  wave- wave  interactions  and  topographic  form- 
drag.  Their  model  is  extended  to  include  a  zonally  asymmetric  vorticity  source,  which 
is  shown  to  interact  with  the  perturbation  field  to  produce  zonally  rectified  wave  fluxes 
that  dramatically  alter  the  Hopf  bifurcations  from  stationary  solutions  to  low  frequency 
oscillations.  The  frequency,  intensity,  and  general  character  of  these  oscillations  are  shown 
to  depend  crucially  upon  the  phasing  and  relative  strength  of  the  forcings. 
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Variational  Data  Assimilation  Experiments  with  Variable 
Resolution  Finite-Element  Shallow- Water  Equations  Model 

by 

Keyun  Zhu 

Supercomputer  Computations  Research  Institute 
Florida  State  University,  Tallahassee,  FL  32306-4052 

I.  M.  Navon 

Department  of  Mathematics,  Geophysical  Fluid  Dynamics  Insitute  and 

Supercomputer  Computations  Research  Institute 
Florida  State  University,  Tallahassee,  FL  32306-1307 

and 

X.  Zou 

Supercomputer  Computations  Research  Institute 
Florida  State  University,  Tallahassee,  FL  32306-4052 

May  28,  1998 

Submitted  to  Monthly  Weather  Review 
Abstract 

The  adjoint  model  of  a  finite-element  shallow-water  equations  model  was  obtained 
with  a  view  to  calculate  the  gradient  of  a  cost  functional  in  the  framework  of  using  this 
model  to  carry  out  variational  data  assimilation  (VDA)  experiments  using  optimal  control 
of  partial  differential  equations. 

The  finite-element  model  employs  a  triangular  finite-element  Galerkin  scheme  (Navon, 
1979).  The  derivation  of  the  adjoint  of  this  finite-element  model  involves  overcoming 
specific  computational  problems  related  to  obtaining  the  adjoint  of  iterative  procedures 
for  solving  systems  of  non-symmetric  linear  equations  arising  from  the  finite-element  dis¬ 
cretization  and  dealing  with  irregularly  ordered  discrete  variables  at  each  time-step. 

The  correctness  of  the  adjoint  model  was  verified  at  the  subroutine  level,  and  was 
followed  by  a  gradient  check  conducted  once  the  full  adjoint  model  was  assembled.  VDA 
experiments  were  performed  using  model- generated  observations.  In  our  experiments, 
assimilation  was  carried  out  assuming  that  observations  consisting  of  a  full-model-state 
vector  are  available  at  every  time  step  in  the  window  of  assimilation.  Successful  retrieval 
was  obtained  using  the  initial  conditions  as  control  variables,  involving  the  minimization 
of  a  cost  function  consisting  of  the  weighted  sum  of  difference  between  model  solution  and 
model  generated  observations. 

An  additional  set  of  experiments  was  carried  out  aiming  at  evaluating  the  impact 
of  carrying  out  VDA  involving  variable  mesh  resolution  in  the  finite-element  model  over 
the  entire  assimilation  period.  Several  conclusions  are  drawn  related  to  the  efficiency  of 
VDA  with  finite  element  discretization  with  variable  horizontal  mesh  resolutions  and  the 
transfer  of  information  between  coarse  and  fine  meshes. 


The  Global  Atmospheric  Response  to  Localized 
Low  Frequency  Forcing:  Zonally  Averaged  Basic  States 

by 

Long  Li  and  Terrence  R.  Nathan 
Atmospheric  Science  Program 
Department  of  Land,  Air  and  Water  Resources 
University  of  California,  Davis,  CA  95616 
submitted  to 

Journal  of  Atmospheric  Sciences 
August  1993 

Abstract 

The  remote  response  to  localized,  low  frequency  forcing  is  examined  using  a  linearized, 
non-divergent  barotropic  model  on  a  sphere.  Zonal  mean  basic  states  characterized  by 
solid  body  rotation,  critical  latitudes,  reflective  surfaces,  or  wave  guides  are  considered. 
For  slowly  varying  basic  states  in  which  the  forcing  is  located  entirely  within  the  tropical 
easterlies,  a  WKB  and  ray  tracing  analysis  shows  that  in  sharp  contrast  to  stationary 
waves,  low  frequency  waves  can  propagate  through  the  tropical  easterlies  into  middle  lat¬ 
itudes.  This  propagation  is  most  pronounced  for  wave  periods  of  10-20  days.  Based  on 
the  principle  of  conservation  of  wave  action,  it  is  shown  that  as  the  meridional  scale  of  the 
wave  decreases  (increases)  the  wave  amplitude  must  increase  (decrease).  An  eigenanalysis 
shows  that  although  both  eastward  and  westward-propagating  normal  modes  are  gener¬ 
ated  by  the  localized,  low  frequency  forcing,  the  forcing  projects  more  strongly  onto  the 
westward-propagating  modes,  which  are  favored  for  meridional  propagation.  These  results 
strongly  suggest  that  westward-propagating,  low  frequency  disturbances  generated  within 
the  tropics  may  have  a  more  important  impact  on  long-term  weather  prediction  in  middle 
latitudes  than  previously  thought. 
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III.  SUMMARY  OF  ACCOMPLISHMENTS 
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1.  Low  Frequency  Dynamics  and  Predictability 


Topographic  (i.e.,  form-drag)  instability  is  one  of  the  mechanisms  proposed  as  playing 
an  important  role  in  atmospheric  low  frequency  variability  (e.g.,  Frederiksen  1983;  Legras 
and  Ghil  1985;  Nathan  1988, 1989;  Howell  and  Nathan  1990;  Frederiksen  and  Bell  1987;  Jin 
and  Ghil  1990;  Tribbia  and  Ghil  1990;  O’Brien  and  Branscome  1991).  Drs.  Barcilon  and 
Nathan  conjecture  that  understanding  block  formation  and  demise  may  be  linked  to  un¬ 
derstanding  low  frequency  dynamics.  As  originally  shown  by  Charney  and  DeVore  (1979), 
alteration  of  the  zonal  mean  momentum  will  likely  play  an  important  role  in  topographi¬ 
cally  induced  instabilities.  This  point  is  underscored  by  Nathan  (1992),  who  analytically 
demonstrated  within  the  context  of  the  two-layer  baroclinic  model  that  steady,  zonally 
varying  potential  vorticity  (PV)  sources  can  interact  with  the  wave  field  to  alter  the  zonal 
mean  flow,  resulting  in  dramatic  changes  in  the  linear  topographic  instability  properties 
of  the  system.  However,  at  finite-amplitude  only  steady  solutions  were  obtained;  the  ab¬ 
sence  of  perpetual,  low  frequency  oscillations  was  attributed  to  the  exclusion  of  wave- wave 
interactions  within  the  model. 


Under  AFOSR  Grant  49620-92- J-0498,  Drs.  Barcilon  and  Nathan  are  studying  the 
impact  of  steady,  zonally  varying  vorticity  forcing  on  topographic  instability  and  finite 
amplitude,  low  frequency  oscillations.  The  model  consists  of  a  viscous,  barotropic  fluid  on 
a  zonally  periodic,  midlatitude  /3-plane  channel,  bounded  above  by  a  flat,  rigid  boundary, 
below  by  a  bottom  topography,  and  laterally  by  sidewalls  at  y  =  0,  ir.  For  a  steady, 
meridionally  sheared  basic  flow  U(y )  the  nondimensional  equations  governing  the_quasi- 
geostrophic  dynamics  of  a  superimposed  disturbance,  ip(x,y,t)  and  mean  field,  rp(y,t), 
are: 


(jg  +  V-i^)  VV  +  ifi  -  Vn)j£  +  V2*)  +  j(<p,  hB) 


=  -7V2[^>  -  »(»)]  -  V 


Tr  dhB 


dx 


d  d 2ip 

dtlw 


=  _  i.  (WhA  _  ±  +  7 

1  dy*  dy\dxB)  dy\dx  W )  7 


(la) 

+  qB  , 

,\  a2* 

(16) 

with  ifix  =  ipt  =  0  at  y  =  0,  7T.  In  these  equations  hs{x,y)  is  the  bottom  topography, 
9B(x,y)  is  a  steady,  zonally  varying  vorticity  source,  'k(y)  is  an  imposed  perturbation 
zonal  momentum  drive  and  7  measures  the  Ekman  pumping  strength  on  the  upper  and 
lower  boundaries. 


Equation  (16)  indicates  that  changes  in  the  mean  flow  vorticity  are  due  to  damping, 
topographic  form  drag,  convergence  of  zonally  rectified  vorticity  fluxes  and  externally 
imposed  drive. 

Nathan  and  Barcilon  (1993)  let  U(y)  =  Ug(y),  where  U  is  a  characteristic  zonal  veloc¬ 
ity?  hB(x,y)  =  H(y){helkx  +  *}  be  the  topographic  forcing  and  </b(x ,  y)  =  Q{y){qtlkx  +  *} 
be  the  vorticity  forcing.  The  forcing  amplitudes  h  and  q  are  assumed  complex,  thus  allowing 
for  a  phase  shift  between  the  forcings.  Then  they  let  U  =  where  |y|  1  measures 
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the  departure  from  linear  resonance  produced  by  the  perturbation  zonal  momentum  source 
’i'(y)  and  introduce  the  long  time  scale  r  =  gt;  the  dissipation  and  forcing  terms  are  ordered 
as  follows  7  =  e2R;  ha  =  c3h(x,y);  qa(x,y)  =  e3g(x,y);  (V»)  =  s2\P<2>(y);  e2  =  0(g). 
They  write 

iJ>(x,y,T)  =  etfi(y){A{r)etkx  +  *}  +  e2  F(y){A2  (r)ei2kx  +  *}  +  eVo2>(y>r)  •  (2) 

For  arbitrary  phasing  between  the  topographic  and  vorticity  forcings,  the  governing 
equations  are  non-separable  in  y  and  r.  To  simplify  the  analysis,  Nathan  and  Barcilon 
assume  that  the  vorticity  forcing  is  maximized  either  90°  upstream  or  90°  downstream  of 
the  topography;  omitting  the  details  reported  in  their  paper,  the  governing  equations  are 
now  of  the  form 


— — |-  RA  +  iSA  +  iu;|.A|2  +  iaZA  —  ibhr  +  =  0  ,  (3a) 

ar 

and 

^  +  RZ  =  -A,  ,  (36) 

dr 

where  all  the  coefficients  are  real  and  defined  in  their  paper. 

The  second  term  in  (3a)  represents  linear  (Rayleigh)  damping  and  the  third  term 
represents  the  linear  frequency  correction  arising  from  the  0(g)  departure  of  the  basic 
flow  from  resonance.  The  first  three  tenns  alone  yield  a  damped,  slowly  propagating 
Rossby  wave.  The  fourth  term  is  the  nonlinear  frequency  correction  arising  from  wave- 
wave  interaction  and  is  essential  to  the  generation  of  low  frequency  oscillations.  The  fifth 
term  represents  wave-mean  flow  interaction.  The  last  two  terms  are  the  forcings,  which 
may  augment  or  oppose  each  other  depending  on  the  spatial  structures,  strengths,  and 
phasing.  The  strength  of  the  topographic  forcing,  measured  by  6,  is  proportional  to  Ui  k 
implying  that  the  topographic  forcing  dominates  over  the  vorticity  forcing  when  the  zonal 
winds  are  strong  and  the  zonal  wave  scale  short. 

To  illustrate  the  role  of  vorticity  forcing  on  the  topographically  induced  Hopf  bifur¬ 
cation,  eqs.  (3a)  and  (36)  are  written  as 


Xi  +  RXi  -  (6  +  X3)X2  -  w(X2  +  X2)X2  =  0  ,  (4a) 

X  +  RX 2  —  (6  +  X3)X\  -+•  &(Xi  +  Xi)  =  1  -f  gi  ,  (4b) 

X  -I-  RX3  =  —cX2  (4c) 

where  X\  +  iX2  —  A/bhr,  X3  =  aZ,  c  =  abh 2,  u>  =  u)(bhr)3,  and  g,  =  —b2qi/bh. 


The  dynamical  system  described  by  (4)  was  integrated  numerically  using  a  fourth  order 
Runge-Kutta  scheme;  a  nondimensional  time  step  of  0.05  generally  proved  sufficient.  For 
the  cases  presented  below,  Nathan  and  Barcilon  have  varied  6  and  qi  and  set  the  remaining 
parameters  and  initial  conditions  equal  to  those  used  by  Jin  and  Ghil  (1990)  in  their  Figs. 
4-8:  c  =  1.0,  u  =  0.5,  R  =  0.1,  ^(0)  =  A'2(0)  and  X3(0)  =  -4.8. 
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Figure  1  summarizes  the  types  of  solutions  that  emerge  from  (4)  for  different  combi¬ 
nations  of  6  and  qi.  The  steady  solutions  are  denoted  by  5.  The  oscillatory  solutions  are 
generally  of  two  types.  The  first  type,  denoted  by  P,  has  been  identified  by  Jin  and  Ghil 
(1990)  and  is  characterized  by  a  single,  low  frequency  oscillation,  which  falls  within  the 
intraseasonal  time  scale.  The  second  type,  denoted  by  T,  only  arises  for  qi  >  0,  i.e.,  when 
the  vorticity  forcing  augments  the  topographic  forcing.  The  T-solutions  are  characterized 
by  a  very  long  oscillation  within  which  is  embedded  two  distinctly  different  regimes  of  be¬ 
havior:  a  damped,  low  frequency  amplitude  oscillation  in  which  the  wave  field  propagates 
around  the  domain,  and  a  state  where  the  wave  field  remains  quasi-steady.  Only  steady 
solutions  exist  for  |6j  sufficiently  large. 


Figure  1:  Mapping  of  the  steady  (5),  low  frequency  oscillation  (P),  and  two  regime 
oscillation  (T)  solutions  that  arise  for  various  combinations  of  6  and  qi,  where  8  is 
proportional  to  the  departure  from  topographic  resonance  and  q,  is  proportional  to  the 
amplitude  of  the  vorticity  forcing.  The  solid  lines  separates  the  steady  and  oscillatory 
solutions  and  the  dashed  line  separates  the  single  and  two  regime  solutions,  pw 


Topographically  induced  low  frequency  oscillations  in  the  presence  of  a  steady,  zonally 
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varying  vorticity  source  have  been  studied.  The  source  has  been  shown  to  interact  with  the 
perturbation  field  to  produce  zonally  rectified  wave  fluxes  that  change  the  zonal  momentum 
balance  of  the  system.  This  change  provides  a  mechanism  for  coupling  the  topographic  and 
vorticity  forcings,  leading  to  a  dramatic  alteration  in  the  Hopf  bifurcation  from  stationary 
to  low  frequency,  finite  amplitude  oscillations.  In  particular,  Nathan  and  Barcilon  have 
shown  that  the  vorticity  forcing  is  capable  of  generating  low  frequency  oscillations  that 
remain  strictly  upstream  of  the  topographic  ridge.  For  sufficiently  large  vorticity  forcing, 
the  Hopf  bifurcation  is  destroyed,  irrespective  of  the  relative  phase  of  the  forcings.  Nathan 
and  Barcilon  have  obtained  solutions  that  exhibit  two  distinctly  different  types  of  regimes. 
One  regime  is  a  westward-propagating,  low  frequency  oscillation  that  has  the  character  of 
a  free  wave.  The  second  regime  is  characterized  by  a  quasi-steady  state  in  which  the  wave 
ridge  remains  upstream  of  the  topographic  ridge.  Details  are  contained  in  the  paper  by 
Nathan  and  Barcilon  (1993). 


2.  Use  of  Empirical  Orthogonal  Functions  (EOFs)  as  Basis  Functions 


a.  Introduction 

Lorenz  (1956)  first  used  empirical  orthogonal  functions  (EOFs)  as  basis  functions  in  the 
analysis  of  atmospheric  data  in  place  of  tabulated  functions,  such  as  spherical  harmonics. 
Typically,  the  number  of  EOFs  required  to  represent  95%  of  the  variance  of  a  meteorolog¬ 
ical  variable  is  about  one  third  the  number  of  degrees  of  freedom  found  in  conventional 
numerical  weather  prediction  models. 

One  of  the  disadvantages  of  EOFs  is  their  lack  of  well  defined  functional  forms  and 
the  associated  difficulty  in  dealing  with  these  functions  in  mathematical  equations.  Schu¬ 
bert  (1985)  showed  that,  with  the  use  of  spherical  harmonics,  it  is  possible,  in  part,  to 
overcome  this  limitation  by  formulating  the  EOFs  in  the  spectral  domain  and  using  spher¬ 
ical  harmonic  identities  to  simplify  the  mathematical  expressions.  Schubert  developed  a 
barotropic  EOF  forecast  model  that  conserves  either  energy  or  enstrophy  for  the  truncated 
EOF  set;  both  energy  and  enstrophy  are  conserved  if  the  full  set  of  EOFs  is  used.  Using  40 
EOFs  with  an  equivalent  barotropic  model,  he  found  equally  good  correlations  between  the 
observed  and  EOF-modeled  stream  function  tendencies  as  were  found  between  observed 
and  spectrally-modeled  tendencies  using  68  spectral  coefficients. 

Graduate  student  Y.  Chang,  in  collaboration  with  Drs.  Barcilon  and  Pfeffer,  is  applying 
the  EOF  methodology  to  a  realistic,  but  idealized,  two-layer,  baroclinic  model  to  under¬ 
stand  the  advantages  and  drawbacks  of  EOF  modelling  in  numerical  weather  prediction. 
The  EOF  is  based  on  very  long  model  data  sets  generated  using  Dr.  Whitaker’s  spectral 
model  developed  while  he  was  a  graduate  student  in  our  program  (Whitaker  1990). 
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b.  Summary  of  research  in  numerical  prediction  using  EOFs 

The  words  italicized  below  correspond  to  components  in  the  block  diagram  displayed 
in  Figure  2. 

The  two-layer  baroclinic  hemispherical  primitive-equation  (spectral  transform)  model 
has  been  used  to  generate  100  years  of  data.  After  subtracting  the  100  year  time  mean  from 
the  data  set,  Mr.  Chang  computed  a  set  of  real  EOFs  for  each  field  in  each  layer,  which 
was  used  to  construct  the  EOF  model  (Basis  Functions  and  EOF  Model).  To  examine  the 
significance  of  the  data  set,  he  used  the  dominant  variance  rule  (Rule  N)  and  performed  the 
Random  Test.  He  also  examined  Eigenvalue  Resolution  which  provides  an  approximation 
to  the  sampling  error  of  the  eigenvalues  and  eigenvectors  of  the  covariance  matrix.  From 
the  Random  Test  and  Eigenvalue  Resolution  he  obtained  a  measure  of  the  number  of  EOFs 
needed. 


Figure  2. 
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The  following  presents  in  greater  detail  the  work  which  is  summarized  in  the  preceding 
paragraph. 


c.  Data  and  data  analysis 

c.l  Data:  Under  the  supervision  of  Drs.  Barcilon  and  Pfeffer,  Mr.  Chang  analyzed  100 
years  of  data  for  stream  function,  velocity  potential  and  the  thickness  for  each  layer.  A 
spherical  harmonic  (K^”)  expansion  was  then  performed  with  T31  truncation.  The  data 
consist  of  twice  daily  values  of  three  variables  for  each  layer  in  one  hemisphere  of  the 
model.  The  bottom  topography  is  a  zonal  wave  number  two  mountain  with  peaks  at  45 
degrees  latitude  and  the  mountain  height  vanishing  at  the  north  pole  and  at  the  equator; 
its  amplitude  is  2  km.  Figure  3  shows  upper  layer  mean  stream  function  and  the  variance 
of  the  meridional  velocity  v  for  100  years  of  data. 


c.2  Empirical  Orthogonal  Functions:  Chang,  Barcilon  and  Pfeffer  expand  the  data  in 
spherical  harmonics  and  write 

.F(A,  /x,  t)  =  ^  ]  F-y(t)Y^( A,  /x)  (3.2.1) 

i 

where  F7  sure  the  (complex)  spectral  transforms  of  amplitude  (Platzman).  We  then  seek 
an  orthogonal  function  expansion  of  Fy(t)  of  the  form  (Schubert,  1985) 

*V(0 -£*»(*)*( 7),  (3.2.2) 

k 

Ekiy)  are  the  complex  EOFs  in  the  wave  number  domain  which  form  a  complete  set,  and 
are  found  by  solving  the  appropriate  eigenvalue/vector  problem  for  the  spectral  coefficient 
covariances.  The  principal  components  (PCs),  r*(f),  are  given  by 

Zk(t)  =  '£tEi(  7)F7.  (3.2.3) 

7 

The  upper  layer  streamfunction  EOF  structures  are  shown  in  Figure  4.  EOFl  has  a 
dominant  zonal  mean  and  wave  number  2  pattern,  and  accounts  for  up  to  11.4%  of  the 
total  variance.  It  contains  many  of  the  features  of  the  vn  plots  (Figure  3).  EOF2  and 
EOF3  are  associated  with  a  traveling  wave  #5,  and  account  for  15.4%  of  the  variance. 
EOF4  and  EOF5  produce  a  traveling  wave  #6  and  represent  11.1%  of  the  variance.  The 
NS  oriented  dipole  structures  over  the  mountain-valley  regions  in  EOF6  are  believed  to 
be  related  to  blocking.  EOF7  has  important  contributions  from  zonal  wave  number  2  in 
lower  latitudes,  zoned  wave  number  4  in  midlatitudes  and  zonal  wave  number  0  in  higher 
latitudes.  EOF8  is  made  up  predominantly  of  zonal  wave  number  1.  This  is  reflected  in 
the  spatial  pattern  which  is  composed  of  2  large  anomaly  centers  of  opposite  sign  on  the 
two  mountain  ridges. 
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LATER  2  PS!  MEAN 


Fig.  3  a.  Upper  layer  mean  stream  function,  b.  Variance  vn 
of  meridional  velocity.  In  the  figures,  zero  deg  longitude  is  at  3  o’clock 
and  the  mountain  peaks  are  at  45°  and  135°  longitude. 
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Fig.  4  EOFs  of  streamfunction  $  for  the  upper  layer 


c.S  Selection  Rule  for  Principal  Component  Analysis:  A  given  field  Z  =  z(x,t)  can  be 
represented  as 

P 

*(*»*)  ’ 
j=i 


(3.3.1) 


where  2  =  1, P\t  —  1, N,  and  the  sequence  of  eigenvalues  is  such  that  Ai  >  ...A  p  >  0. 
Generally,  after  a  certain  index  P\  the  magnitudes  of  the  A j  drop  abruptly  and  remain 
relatively  small;  therefore,  only  the  first  P'  eigenvalues  may  be  considered  important  in 
the  representation  of  the  total  variance  of  the  data  set  (3.1.1).  The  remaining  eigenvalues 
contain  negligible  variance  and  are  difficult  to  differentiate  from  those  generated  by  noise. 


Chang,  Bardlon  and  Pfeffer  used  data  sets  of  length  100  years  to  test  the  significance 
of  the  truncated  data  set  (Preisendorfer  et  al.,  1981).  They  used  the  asymptotic  behavior 
of  random  reference  data  test  and  compared  the  noise  eigenvalues  to  the  data  eigenvalues. 
The  95%  points  on  the  cumulative  distribution  of  the  noise  process  (small  dots)  are  shown 
in  Figure  5.  The  data  eigencurve  (large  dots)  is  shown  cutting  through  the  noise  eigencurve. 
The  data  eigenvalues  are  no  longer  considered  significantly  greater  than  noise  eigenvalues 
when  found  beneath  that  threshold.  For  the  stream  function  field,  the  first  50  PCs  have 
variability  significantly  larger  than  would  be  expected  from  a  noise  process.  These  modes, 
when  added,  explain  about  90%  of  the  variance. 


d.  The  EOF  model 


d.l  Model  formulation:  Each  quantity  is  expressed  in  terms  of  spherical  harmonics  and 
the  orthonormality  of  Yy  is  used  to  isolate  the  tendency  equation  for  each  component.  The 
governing  equations  become 

Ct<7(0  -  W,  X,  Z) , 

<hX'n,i(t)  «  <?(*  X,  Z)  , 

Z»,y(t)  -  *(*,  X,  Z)  , 

where  rj>  is  the  streamfunction,  x  is  the  velocity  potential  and  Z  is  the  layer  thickness.  The 
dot  represents  a  time  derivative.  The  spectral  quantities  are  expanded  in  terms  of  EOFs 
as 


k 

x»,«(o = y.  4 ”  ww . 

k 

zu*)  =  E  Ckn)W?z'n(6) , 


where  Aj^(<)  and  C*n*(t)  are  the  real  PCs. 
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Fig.  5  Eigenvalues  for  EOFs  1  to  100,  and  95%  confidence  limits. 


After  some  manipulation,  we  obtain  for,  say  the  vorticity  equation,  an  equation  of  the 
form: 

Ajn)(t)  =  NL  4-  Rotation  +  Mean/EOF  +  Dissipation, 

where  “Rotation”  denotes  a  collection  of  terms  involving  the  Coriolis  and  0  parameters, 
“Mean /EOF”  denotes  terms  associated  with  the  mean  flow  and  EOF  interactions  and  the 
other  remaining  symbols  have  standard  meanings. 

d.t  Variance  of  the  EOF  model  terms:  The  standard  deviations  of  the  various  terms  in 
the  model  are  calculated  as  a  function  of  the  PCs  using  a  one  year  data  set.  The  Spectral 
Transform  Model  (ST  Model)  tendency  is  computed  from  a  centered  finite  difference  for¬ 
mula.  All  other  terms  on  the  right  hand  side  of  that  equation  are  computed  by  inserting 
the  sampling  values  of  the  PCs  into  the  EOF  model  terms. 

Figure  6  shows  the  comparison  between  the  standard  deviation  of  the  EOF  Model 
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tendency  and  that  of  the  ST  Model  for  the  vorticity  equation  for  upper  layer.  Th'ire 
is  a  very  good  agreement  for  the  PCs  of  the  50  EOF  Model.  In  the  upper  layer,  the 
dominant  terms  are  the  mean  flow/EOF  interactions  and  the  Rotation  term.  These  two 
terms  are  much  larger  than  the  tendency  term  in  some  large  zonal  wave  number  PCs,  but 
are  highly  negatively  correlated  so  that  their  sum  has  a  much  smaller  standard  deviation. 
Two  large  peaks  occur  in  the  mean  flow/EOF  term,  one  around  the  2nd  ,  3rd  ,  4th  and  5th 
PCs  and  others  at  the  10th  and  11th  PCs.  The  spatial  distribution  of  the  corresponding 
EOF  patterns  reveals  very  strong  traveling  waves  with  zonal  wave  number  5-7.  The 
EOF  interaction  terms  are  smaller  than  the  tendency  for  all  components  with  the  largest 
contribution  coming  from  the  9th  component.  The  vertical  advection  term  is  at  least  one 
order  of  magnitude  smaller  than  the  tendency  term  for  all  the  PCs. 

STANDARD  DEVIATION  OF  MODEL  TERMS 


d.S  Tendency  Correlations:  A  simple  measurement  of  the  forecast  skill  of  the  EOF 
Model  is  to  compute  the  correlations  between  the  ST  Model  tendency  and  the  EOF  Model 
tendencies  for  each  PC.  The  method  we  used  here  is  identical  to  that  of  the  previous 
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section. 

Figure  7(a)  shows  the  correlations  as  a  function  of  PC  for  the  models  truncated  at  50 
PCs.  The  correlations  are  presented  between  ST  Model  tendency  and  Mean/EOF,  EOF 
interactions  and  Rotation  terms  of  the  EOF  Model.  In  the  upper  layer,  one  finds  the 
Rotation  term  and  Mean/EOF  term  are  negatively  correlated  for  most  PCs.  They  have 
almost  the  same  magnitude  but  opposite  sign.  The  EOF  interaction  term  is  important  in 
the  upper  layer.  It  has  a  correlation  of  0.5  for  the  2nd  and  3rd  PCs.  This  high  correlation 
shows  that  the  main  contribution  to  the  tendency  comes  from  the  nonlinear  interaction 
terms. 

Figure  7(b)  shows  the  complete  50  PCs  EOF  Model  tendency.  In  the  upper  layer, 
seven  PCs  have  correlations  greater  than  0.8;  these  are  the  2nd  ,  3rd  ,  5th  ,  6th  ,  8th  ,  9th  and 
10th .  There  is  a  big  drop  in  correlation  amplitude  when  the  EOF  interaction  terms  are 
neglected.  The  results  show  a  trend  toward  lower  correlations  with  increasing  PC  number. 


e.  Future  Work:  See  Figure  1  for  a  sketch  of  the  future  work  which  we  propose  to  undertake 
next  year. 


3.  Laboratory  Model  of  Moist  Convection  and  Parameterization  in  NWP 

Dr.  Ruby  Krishnamurti  is  completing  a  paper  entitled  “Generation  of  low  frequency 
modes  in  turbulent  Rayleigh-Bemard  convection”,  which  will  be  submitted  shortly  to  J. 
Fluid  Mech..  This  paper  presents  a  mathematical  model  that  accounts  for  a  major  new 
experimental  result,  namely  the  discovery  that  steady  forced  turbulent  convection  cam 
spontaneously  generate  traveling  organized  structures,  giving  rise  to  a  long  period  oscilla¬ 
tion  superimposed  on  the  high  frequencies  usually  associated  with  turbulent  convection. 

Dr.  Krishnamurti  has  also  designed  and  constructed  an  apparatus  for  measuring  heat 
and  momentum  fluxes  in  an  annulus  of  convecting  fluid.  In  particular,  the  fluid  is  stably 
stratified  by  heating  from  above  and  cooling  from  below,  but  convects  because  of  internal 
heat  sources.  The  location  of  this  internal  heating  is  governed  by  the  fluid  flow  itself: 
the  working  fluid  is  water  containing  a  small  amount  of  thymol  blue  which  is  normally 
used  as  a  pH  indicator.  In  its  neutral  state  it  makes  the  water  appear  orange,  but  with 
an  excess  of  hydroxyl  ions,  it  makes  the  water  dark  blue.  In  the  present  arrangement, 
the  working  fluid  acts  as  a  very  weak  electrolyte,  its  lower  boundary  acts  as  the  positive 
electrode  connected  to  a  9  volt  battery.  Near  this  lower  boundary  then,  the  water  turns 
dark  blue.  Orange  light  from  a  sodium  vapor  lamp  passes  through  the  orange  water  with 
negligible  absorption,  but  is  strongly  absorbed  by  the  blue  fluid  which  is  thereby  heated. 
(The  instability  of  this  heated  layer  is  the  subject  of  a  paper  entitled  “Stability  theory  of 
an  internally  heated  fluid  layer”  by  R.  Krishnamurti  and  S.  Ranganathan.)  Instability  in 
this  heated  layer  causes  blue  fluid  to  rise.  However,  as  long  as  it  remains  blue,  it  continues 
to  absorb  radiation,  continues  to  be  heated  and  continues  to  rise.  Thus  the  internal  heat 
source  is  only  in  the  rising,  (not  in  the  sinking)  plumes. 
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Fig.  7  Correlation  coefficients  as  a  function  of  principal  component, 
a.  Correlations  between  “observed"  tendency  and  each  term  in  the  tendency  equation, 
b.  Correlations  between  “observed”  tendency  and  model  tendency  with  and  without  nonlinear  NL  terms. 


The  bottom  boundary  is  driven  by  a  synchronous  motor,  providing  a  Couette  flow 
with  vertical  shear.  In  this  environment,  the  momentum  transport,  sis  well  as  the  heat 
transport,  by  the  convecting  fluid  will  be  measured.  The  apparatus  should  be  operational 
by  Fall  1993. 


D.  Finite  Element  Model  of  the  Shallow- Water  Equations 


A  copy  of  the  paper  by  Zhu,  Navon  and  Zou  submitted  to  Mon.  Wea.  Rev.  is  being 
sent  under  separate  cover. 
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